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SUMMARY 

By using the theory and numerical procedure developed by the authors in earlier 
publications, the influence of lubricant starvation upon minimtim film thickness in 
starved elliptical elastohydrodynamic conjunctions for low -elastic -modiilus materials 
has been investigated. Lubricant starvation was studied simply by moving the inlet 
boundary closer to the center of the conjunction. The results show that the location of 
the dimensionless inlet boundary m* between the fully flooded and starved conditions 
can be expressed simply as m* = 1 + 1,07 [(R /b)^ H . where R is the 

effective radius of curvature in the rolling direction, b is the semiminor axis of the 
contact ellipse, and p, is the dimensionless minimum film thickness for the fully 

flooded condition. That is, for a dimensionless inlet distance m less than m*, star- 
vation occurs; and for m s m*, a fully flooded condition exists. Furthermore, it has 
been possible to express the minimum film thickness for a starved condition as 

Contour plots of the pressure and film thickness in and around the contact are pre- 
sented for both the fully flooded and starved lubrication conditions. It is evident from 
the contour plots that the inlet pressure contours become less circular and that the film 
thickness decreases substantially as the severity of starvation increases. 

The results presented in this report, when combined with the findings previously re- 
ported, enable the essential features of starved, elliptical, elastohydrodynamic conjunc- 
tions for materials of low elastic modulus to be ascertained. 

*Presented at Joint Lubrication Conference sponsored by the American Society of 
Lubrication Engineers and the American Society of Mechanical Engineers, Minneapolis, 
Minnesota, October 24-26, 1978. 

'Professor of Mechanical Engineering, University of Leeds, Leeds, England. 



INTRODUCTION 


Only in recent years has the complete numerical solution of the isothermal elasto- 
hydrodynamic lubrication (EHL) of elliptical contacts successfully emerged. The anal- 
ysis requires the simultaneous solution of the elasticity and Reynolds equations. The 
authors' approach to the theoretical solution has been presented in two previous publica- 
tions (refs. 1 and 2). The first of these publications (ref, 1) represents an elasticity 
model in which the conjimction is divided into equal rectangular regions with a uniform 
pressure applied over each region. The second (ref. 2) gives a complete approach to the 
solution of the elastohydrodynamic lubrication problem for elliptical contacts. 

The most important practical aspect of the EHL elliptical -contact theory (ref. 2) is 
the determination of the minimum film thickness within the contact. That is, the pre- 
diction of a film of adequate thickness is extremely important for the successful opera- 
tion of machine elements in which these thin, continuous, fluid films occur. In refer- 
ence 3 the fully flooded results obtained from the theory given in references 1 and 2 are 
presented. A fully flooded condition is said to exist when the inlet distance of the con- 
junction ceases to influence in any significant way the minimum film thickness. In ref- 
erence 3 the influence of the ellipticity parameter and the dimensionless speed, load, 
and material parameters on minimum film thickness was investigated. Thirty -four 
cases were used in obtaining the fully flooded minimum -film -thickness formula. 

In reference 4 the basic theory developed in references 1 and 2 was used to study the 
effect of lubricant starvation on the pressure and film thickness within the conjunction. 

A simple expression for the dimensionless inlet bovmdary distance was obtained. This 
inlet boundary distance defines whether a fully flooded or a starved condition exists in 
the contact. Fifteen cases, in addition to three presented in reference 3, were used to 
obtain simple expressions for the minimum and central film thicknesses in a starved 
conjunction. 

The work presented in references 1 to 4 is for a material of high elastic modulus 
(e, g. , steel). The work presented in reference 5 and in the present paper is for a mate- 
rial of low elastic modulus (e. g. , nitrile rubber). For such materials the distortions 
are large even with light loads. Another feature of the EHL of low -elastic -modulus ma- 
terials is the negligible effect of pressure on the viscosity on the lubricating fluid. 
Engineering applications in which elastohydrodynamic lubrication is important for low- 
elastic-modiilus materials include seals, human joints, tires, and elastomeric -material 
machine elements. 

The problem of fully flooded line contacts has been solved theoretically for low- 
elastic -modulus materials by Herrebrugh (ref. 6), Dowson and Swales (ref. 7), and 
B^lin and Archard (ref, 8), The solutions of references 6 and 7 were obtained numer- 
ically and are based on simxiltaneous solutions of the hydrodynamic and elasticity equa- 
tions; the analytical solution of reference 8 relied on the assumption of a simplified form 
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for the film shape in the contact region. Biswas and Snidle (ref. 9) used the approach of 
reference 8 to solve the point -contact situation. Reference 5 presents, to the best of the 
authors' knowledge, the first complete numerical solution of the problem of fully flooded, 
isothermal, elastohydrodynamic lubrication of elliptical contacts for low -elastic -modulus 
materials. That is, no assumptions were made as to pressure or film thickness within 
the contact, and compressibility and viscous effects were considered. 

The present paper uses the basic theory developed in references 1 and 2 in extendii^ 
the work on materials of low elastic modtilus given in reference 5 by studying the effect 
of lubricant starvation on pressure and film thickness within the conjunction. The ap- 
proach used here is similar to that used in reference 4 in studying starvation effects for 
materials of high elastic modulus. That is, the objective of this work is to provide a 
simple expression for the dimensionless inlet boundary distance. This inlet boundary 
distance defines whether a fully flooded or a starved condition exists in the conjunction. 
An additional objective is to develop simple expressions for the minimum film thickness 
under the starvation condition. Eighteen cases, in addition to three presented in refer- 
ence 5, were used in obtainii^ the starvation results. To show more fully what occurs 
in going from a fully flooded condition to a lubricant starvation condition, contour plots 
of pressure and film thickness in and around the contact are presented. 


SYMBOLS 


a 

b 

E 

E’ 


semimajor axis of contact ellipse 
semiminor axis of contact ellipse 
modulus of elasticity 


1 , 1 - 4 

Eg 


F normal applied load 

G dimensionless material parameter, E*/Pjy g^g “E* 

Hmin dimensionless minimum film thickness, 

Hmin -p dimensionless minimum film thickness for a fully flooded conjunction, 

Vln, f/^x 

<3 dimensionless minimum film thickness for a starved conjunction, h^.„ q/R^ 

xxixxx^ inxn^ id 

h film thickness 
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k 

m 

m' 

P 

P 


P 


IV, as 


R 

r 

U 

u 

W 

xX 
y Y 

a 

% 

V 


ellipticity parameter, a/b 
dimensionless inlet distance 

dimensionless inlet distance at boundary between fully flooded and starved 
conditions 

dimensionless pressure, p/E* 
pressure 

asymptotic isoviscous pressure 
effective radius of ctirvature 
radius of curvature 

dimensionless speed parameter, urj./E'R 

mean surface velocity in x-direction, i (u^ + Ug) 

2 

dimensionless load parameter, F/E'R^ 
coordinate system defined in paper 

pressiire -viscosity coefficient 
atmospheric viscosity 
Poisson's ratio 


Subscripts: 

A solid A 

B solid B 

F fully flooded conjunction 

H materials of high elastic modulus 

min minimum 

S starved conjunction 

X, y coordinate system defined in paper 


ANALYSIS 


Boimdary Between Fully Flooded and Starved Conditions 


The nodal structure used in obtaining the fully flooded results of reference 5 is 
shown in figure 1. In this figure the coordinate X is made dimensionless with respect 


4 



to the semiminor axis b of the contact ellipse, and the coordinate Y is made dimen- 
sionless with respect to the semimajor axis a of the contact ellipse. The ellipticity pa- 
rameter k is defined as the semimajor axis divided by the semiminor axis of the con- 
tact ellipse (k = a/b). The coordinates X and Y are made dimensionless such that the 
Hertzian contact ellipse becomes a Hertzian circle regardless of the ellipticity param- 
eter. This Hertzian contact circle is shown in figure 1 with a radius of unity. The pres- 
sure is assTuned to be ambient at the edges of the computing area. In figure 1 the di- 
mensionless inlet distance m, which is equal to the dimensionless distance from the 
center of the contact to the inlet edge of the computing area, is shown. 

Lubricant starvation can be studied simply by reducing the dimensionless inlet dis- 
tance m. A fully flooded condition is said to exist when the dimensionless inlet distance 
ceases to influence in any significant way the minimum film thickness. When starting 
from a fully flooded condition as shown in figure 1 and decreasing m, the value at which 
the minimum film thickness first starts to change is called the fully flooded - starved 
boundary and is denoted by m* , Therefore, lubricant starvation was studied by using 
the basic elastohydrodynamic lubrication theory developed in reference 2, considering 
solid materials of low elastic modulus, and observing the effect of reducing the dimen- 
sionless inlet distance. 

Table I shows how changing the dimensionless inlet distance affects the dimension- 
less film thickness for three groups of dimensionless load and speed parameters. For 
all the results presented in this paper the dimensionless material parameter G is fixed 
at 0. 4276 and the ellipticity parameter k is fixed at 6. As shown in table I, as the di- 
mensionless inlet distance m decreases and dimensionless minimum film thickness 
Hmin decreases. 

Table H shows how the three groups of dimensionless speed and load parameters af- 
fect the location of the dimensionless inlet boundary m* . Also given in this table are 
the correspondii^ values of the dimensionless minimum film thickness for the fully 
flooded condition as obtained by interpolation of the numerical values. The dimension- 
less inlet boundary distance m* shown in table H was obtained by using the data from 
table I when the following equation was satisfied: 


^min,F “ 

^min, F 


0. 03 


( 1 ) 


The value of 0. 03 is used in equation (1) since it was ascertained that the data in table I 
could only be obtained to an accixracy of 3 percent. 

The general form of the equation that describes how the dimensionless inlet distance 
at the fully flooded - starved botmdary varies is given as 
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( 2 ) 


m 


* 


1 + A 



-.B 


H 


min, F 


The right side of equation (2) is similar to the forms of the equation given in Wolveridge, 
et al, (ref. 10), Wedeven, et al. (ref, 11), and reference 4. By using the data obtained 
from table I, the dimensionless inlet boundary distance for materials of low elastic mod- 
\3lus can be written as 


m* =1 + 1.07 



-,0. 16 


^min,F 


(3) 


p 

The coefficient of determination r for these results was good at 0. 960. The value of 
r^ reflects the fit of the data to the resulting equation: 1 being a perfect fit, and zero 
the worst possible fit. A fully flooded condition exists when m > m* , and a starved 
condition when m < m* , From reference 4 the dimensionless inlet boimdary distance 
for materials of high elastic modulus can be written as 


{m*}jj = 1 + 3.34 



-,0.56 


H 

min,F 


(4) 


Comparing equation (3) with equation (4) indicates that the dimensionless inlet boundary 
for materials of Iw elastic modulus m* is less than {m* }jj, therefore implying that 
m* is closer to the center of the contact. 


Film -Thickness Formula 

After the limiting location of the inlet boundary for the fully flooded condition is 
clearly established (eq. (3)), an equation definii^ the dimensionless minimum film thick 
ness for lubricant starvation conditions can be developed. The relation between the di- 
mensionless minimum film thickness in the starved and fully flooded conditions can be 
expressed in general form as 


^min,S ^ ^ / m - 1 \ ° 
®min,F \m* - 1/ 



This form of the equation is exactly the same as that used in reference 4, which was de- 
termined to be the best form after trying a number of other forms. Table IH shows how 
the ratio of the dimensionless inlet distance parameter to the ftilly flooded - starved 
boundary (m - l)/(m* - 1) affects the ratio of minimum film thickness in the starved and 
fully flooded conditions S^^min F* ^ least -squares power curve fit to the 17 

pairs of data points 


I ^min,s\ / m - 1 \ 
\^min,F/. \m* - 


where i = 1,2, . . . , 17 


was used in obtaining the dimensionless minimum film thickness for a starved condition 
for low -elastic -modulus materials 


®^min,S ^min,F 



( 6 ) 
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The coefficient of determination r for these results was excellent at 0. 990. From 
reference 4 the dimensionless film thickness for a starved condition for high -elastic - 
modulus materials was written as 


^^min,S^jj ^^min,F^jj 


m - 1 


vO.25 


* -t 

vm - 1, 


( 7 ) 


The exponents in equations (6) and (7) are quite similar. Therefore, even though the 
fully flooded - starved boundary for materials of low elastic modulus was much less than 
that for materials of high elastic modulus, the starved film thickness equations are quite 
similar. Figure 2 shows the influence of inlet boundary on minimum -film -thickness ra- 
tio for materials of high and low elastic modulus. This figure shows the close agree- 
ment between equations (6) and (7). 

Therefore, whenever m <m"^, where m is defined by equation (3), a lubricant 
starvation condition exists. When this is true, the dimensionless minimxun film thick- 
ness is expressed by equation (6). If m > m* , where m* is defined by equation (3), 
a fully flooded condition exists. The dimensionless minimum film thickness for a fully 
flooded condition p for materials of low elastic modulus was developed in ref- 

erence 5 as 


H 


min, F 


= 7. 43(1 - 0. 85 e"®* 


( 8 ) 



where the dimensionless speed parameter 


U = 


E’R^ 


the dimensionless load parameter 


W = 



and the dimensionless ellipticity parameter 



(9) 


( 10 ) 


( 11 ) 


Simplified Formulas 

To make it easier to calculate the ellipticity parameter k, the elliptic integral of the 
second kind S , and the semiminor axis of the contact ellipse b in equations (3) and (8), 
the approximate expressions from reference 12 are used and are given here as 


where 


11/3 



= 1 + 


5R 


R 

k = 1. 03 1 

,R. 


vO. 64 


( 12 ) 


(13) 


(14) 



( 15 ) 
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1 


1 


(16) 


_1 

R. 


y ^Ay 


^By 


The approximate expressions (eqs. (12) to (14)) enable these terms to be easily calculated 
within 3 -percent accuracy without resorting to charts or numerical methods. 


CONTOUR PLOTS OF RESULTS 

A fuller explanation of what happens in going from a fully flooded to a lubricant star- 
vation condition for materials of low elastic modulus is shown in figures 3 to 6. In fig- 
ure 3 contour plots of dimensionless pressure (P = p/E') are given for group 3 of table I 
and for dimensionless inlet distances of 1. 967, 1. 333, 1. 167, and 1. 033. hi this figure 
and the remainder of the contour plots to be presented, the symbol + indicates the cen- 
ter of the Hertzian contact. It should be remembered that because of the way the coor- 
dinates are made dimensionless, the actual contact ellipse becomes a Hertzian circle re- 
gardless of the ellipticity parameter. The Hertzian contact circle is shown in each fig- 
ure by asterisks. At the top of each figure the contour label and its corresponding value 
are given. The inlet region is to the left and the exit region is to the right. 

The dimensionless pressure contours were kept constant in all parts of figure 3. In 
figure 3(a), with m = 1. 967, a fully flooded condition exists. The contours are circular, 
extending further into the inlet region than into the exit region. As was pointed out in 
reference 5, the pressure spikes found for high -elastic -modulus materials (ref. 3) were 
not in evidence in these solutions for low -elastic -m odulus materials. The lack of a pres- 
sure spike is due to the lack of viscous effects of the fluid in the contact for low -elastic - 
modulus materials. This in turn is due to considerably less presstire being generated in 
a contact with low -elastic -modulus materials than in a contact with high -elastic -modulus 
materials. 

In figure 3(b), with m = 1. 333, a starved condition exists, and the inlet contoxir is 
slightly less circular than that in figure 3(a). In figure 3(c), with m = 1. 167, the star- 
vation is even more severe and the inlet contour is less circular. In figure 3(d), with 
m = 1. 033, the starvation is the most severe and the inlet contour is even less circular. 

In figure 4 contour plots of dimensionless film thickness (H = h/R^^) are given for 
group 3 of table I and dimensionless inlet distances m of 1. 967, 1, 333, 1. 167, and 
1. 033. These film -thickness results correspond to the pressure results in figure 3. 

The central portion of the film -thickness contours becomes more parallel as starvation 
is increased (going from figs. 4(a) to (d)). The minimum -film -thickness area moves to 
the exit region of the contact as starvation is increased. The values of the film - 
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thickness contonrs for the most starved condition (fig, 4(d)) are much lower than those 
for the fully flooded condition (fig. 4(a)). 

Figure 5 more clearly describes these film -thickness restilts. It shows the var- 
iation of the dimensionless film thickness on the X-axis for fovir dimensionless inlet 
distances. The value of Y is held fixed near the axial center of the contact, and the 
dimensionless parameters U and W are held constant as given in group 3 of table I. 

In this figure the central region becomes flatter as starvation becomes more severe. 

Also, in going from a fully flooded condition to a starved condition the film thickness de- 
creases substantially and the location of the minimum film thickness moves closer to the 
exit region. 

Figure 6 shows the variation of the dimensionless film thickness on the Y -axis for 
four dimensionless inlet distances. The value of X is held fixed near the axial center 
of the contact, and the dimensionless parameters U and W are held constant as given 
in group 3 of table I. In this figure, contrary to what is shown in figure 5 for the fully 
flooded condition, the central region is flat and the slopes of the curves increase as star- 
vation increases. 


CONCLUDING REMARKS 

The theory and numerical procedure developed by the authors in earlier publications 
have been used to investigate the influence of lubricant starvation upon minimum film 
thickness in starved elliptical elastohydrodynamic conjimctions for materials of low 
elastic modulus. Lubrication starvation was studied simply by moving the inlet boundary 
closer to the center of the conjunction. The results show that the location of the dimen- 
sionless inlet boundary m* between the fully flooded and starved conditions can be ex- 
pressed simply as 


m* =1 + 1.07 



H 


min,F 


0. 16 


J 


That is, for a dimensionless itilet distance m less than m* , starvation occurs; and for 
m > m* , a fully flooded condition exists. Furthermore, it has been possible to express 
the minimum film thickness for a starved condition as 


^rain,S “®min,F 



10 



ContoTir plots of the pressxire and film thickness in and around the contact have been 
presented for both the fully flooded and starved lubrication conditions. It is evident from 
the contour plots that the inlet pressure contours become less circular and that the film 
thickness decreases substantially as the severity of starvation increases. 

The results presented in this report, when combined with the findings previously 
reported, enable the essential features of starved, elliptical, elastohydrodynamic con- 
junctions for materials of low elastic modulus to be ascertained. 

Lewis Research Center, 

National Aeronautics and Space Administration 
Cleveland, Ohio, April 21, 1978, 

505-04. 
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TABLE I. - EFFECT OF STARVATION ON DIMENSIONLESS 
MINIMUM FILM THICKNESS FOR THREE GROUPS OF 
DIMENSIONLESS SPEED AND LOAD PARAMETERS 


Dimensionless 
inlet distance, 
m 

Group 

1 

2 

3 

Dimensionless load parameter, W 

0.4405X10“^ 

0. 2202x10"® 

0.4405X10"® 

Dimensionless speed parameter, U 

0.5139X10"^ 

0. 1027X10"’^ 

o.sisgxio""^ 

Dimensionless minimum film thickness, 

1.967 

131.8x10"® 

241. 8X10"® 

584. 7X10"® 

1.833 

131.2X10"® 

238. 6X10"® 

572.0X10"® 

1.667 

129.7X10"® 

230.8X10"® 

543. 1X10"® 

1.500 

125.6X10"® 

217.2x10"® 

503.0X10"® 

1.333 

115.9X10"® 

199. 3x10"® 

444. 9X10"® 

1. 167 

98. 11X10"® 

170.4x10"® 

383. 5X10"® 

1.033 

71.80X10"® 

120. 8x10"® 

272. 3X10"® 
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TABLE n. - EFFECT OF DIMENSIONLESS SPEED AND LOAD 


PARAMETERS ON DIMENSIONLESS INLET DISTANCE 


AT FULLY FLOODED - STARVED BOUNDARY 


Group 

Dimensionless 

speed 

parameter, 

U 

Dimensionless 

load 

parameter, 

W 

Dimensionless 

radius 

parameter. 

Dimensionless 
fully flooded 
minimum film 
thickness, 

®min, F 

Dimensionless 
inlet distance 
at fully flooded - 
starved 

boundary, 

* 

m 

1 

0. 5139x10"® 

0.4405X10"® 

19.41 

127. 8x10"® 

1.661 

2 

. 1027X10"'^ 

.2202X10"® 

24.45 

234.5x10"® 

1.757 

3 

. 5139X10"'^ 

.4405X10"® 

19.41 

567.2x10"® 

1.850 


TABLE m. - EFFECT OF DIMENSIONLESS INLET 
DISTANCE ON DIMENSIONLESS MINIMUM - 


FILM-THICKNESS RATIO 


Group 

Dimensionless 
inlet distance, 
m 

Ratio of minimum 
film thickness for 
starved and 
flooded conditions, 

®min, s/**min, F 

Met boxmdary 
parameter, 

(m - l)/(m* - 1) 

1 

1.661 

1 

1 


1.500 

.9828 

.7564 


1.333 

.9069 

. 5038 


1. 167 

.7677 

.2526 


1.033 

.5618 

.0499 

2 

1.757 

1 

1 


1.667 

.9842 

. 8811 


1.500 

.9262 

.6605 


1.333 

. 8499 

.4399 


1. 167 

.7267 

.2206 


1.033 

.5151 

.0436 

3 

1. 850 

1 

1 


1.667 

.9575 

.7847 


1.500 

.8868 

.5882 


1.333 

.7844 

.3918 


1. 167 

.6761 

.1965 


1.033 

.4801 

.0388 





Ratio of minimum fiim thickness for starved and fuiiy 
fiooded conditions, s^^min F 
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Dimensionless 

pressure. 

P =p/E' 

A 0.0123 

B .0121 

C .0117 

D .0110 

E .0100 

F .0085 

G .0060 

H .0030 



(b) m =1.333. 
Figure 3. - Continued. 








Dimensionless 
film thickness, 
H = h/Rx 
A 0.00040 

B .00048 

C .00051 

D .00055 

E .00060 

F .00066 

G .00073 

H .00080 



Figure 4. - Continued. 
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11x10"^ 



Figure 5. - Variation of dimensioniess fiim thickness on 
X-axis for four dimensioniess iniet distances. The 
value of Y is held fixed near axial center of contact 
and the dimensionless parameters U and W are held 
constant as given in group 3 of table I. 


11x10”^ 



Figure 6. - Variation of dimensionless film thickness on 
Y-axis for four dimensionless iniet distances. The 
value of X is held fixed near the axial center of con- 
tact and the dimensionless parameters U and W are 
held constant as given in group 3 of table I. 
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ing direction, b is the semiminor axis of the contact ellipse, and ^ is the dimensionless 

minimum film thickness for the fully flooded condition. That is, for a dimensionless inlet dis- 
tance m less than m*, starvation occurs; and for m & m*, a fully flooded condition exists. 
Furthermore, it has been possible to express the dimensionless minimum film thickness for a 
starved condition as g = p [(m - l)/(m* - 1)] ‘ . Contour plots of the pressure 

and film thickness in and around the contact are shown for the fully flooded and starved lubrica- 
tion conditions . These contour plots show that the inlet pressure contours become less circular 
and that the film thickness decreases substantially as the severity of starvation increases. 
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